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The method of many-particle quantum hydrodynamics has been recently developed, particularly this 
method has been used for an electrically polarized Bose-Einstein condensate. In this paper, we present 
the development of this method for an electrically polarized three dimensional Fermi gas. We derive 
corresponding dynamical equations: equation polarization and equation of polarization current evolution 
as well as the Euler and continuity equations. We study dispersion dependencies of collective excitations 
in a polarized Fermi gas and consider interference of an equilibrium polarization on dispersion properties. 



I. INTRODUCTION 

Electrically polarized ultracold Fermi gases have been 
in the center of attention [Q]] along with the electrically 
polarized Bose-Einstein condensate JU]- JB)]. It is well- 
known that generalization of the Gross-Pitaevskii equation 
has been used for description of the electrically polarized 
Bose-Einstein condensate what was suggested in Ref.s |@]- 
J3]. An analogous non-linear Schrodinger equation can be 
used for study of electrically polarized Fermi gases. How- 
ever, using of the kinetic equation is more prevailing (see 
for example |@], fioh . fllll "). In simple cases the set of 
quantum hydrodynamics (QHD) equations is equivalent to 
corresponding non-linear Schrodinger equation 11211 . but 
it we are interested in studying of evolution of electrical 
dipole moment direction we need more general equation 
set. Such equations can be derived by means of the QHD 
method, which allows to get equation of polarization evolu- 
tion, as it was demonstrated for electrically polarized Bose- 
Einstein condensate lfl3ll - iflfill . This paper is dedicated to 
both the development of the method of many-particle QHD 
for electrically polarized ultracold Fermi gas and studying 
of dispersion of bulk collective excitations inwhere. 

To be certain we will consider fermions with the spin 
1/2. We interesting in dynamic of electric polarization, 
and it's influence on static and dynamic properties of the 
ultracold fermions. We do not consider the dynamic of 
spin of atoms, for simplicity. However, at consideration of 
nonlinear properties of the electrically polarized ultracold 
fermions the role of a spin could be important. The non- 
linear Schrodinger equation is usually used for quantum 
gases studying, so it is important to discuss limits of valid- 
ity of such approximation. In the case of the electrically 
polarized quantum gases, in the corresponding non-linear 
Schrodinger equation a new term is added. However, this 
approximation includes interaction among parallel dipoles 
and influence of this interaction on particles translational 
motion (see Fig. 1 ), and does not account spatial and tempo- 
ral evolution of dipoles direction pictured on Fig. 2. In Ref. 
1 1311 authors developed method accounted spatially inho- 
mogeneous distribution of dipoles directions and it's tem- 
poral evolution. This method includes the potential of in- 
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FIG. 1. (Color online) The figure describes the picture of ultra- 
cold Fermi particles having electric dipole moment and obeying 
to the generalization of the non-linear Schrodinger equation (the 
analog of the Gross-Pitaevskii equation which used for polarized 
Bose-Einstein condensate description). This figure presents sys- 
tem of particles with parallel electric dipole moment. Arrow by 
the circle imitates particle dipole and it's direction. 



teraction unparallel dipoles, and it also includes additional 
equations determining evolution of dipoles moment of a 
volume unit. It was based on the method of many-particle 
QHD. In our previous papers lfl3ll - lfl~7ll we developed 
the method of the many-particle QHD for the electrically 
polarized BEC. This approximation allows to describe as 
the spatial and the temporal evolution of the electric dipole 
moments including evolution of it's directions. Evolution 
of dipoles direction gives influence on translational motion 
of particles due to tensor nature of the dipole-dipole inter- 
action. Equations of the QHD are derived directly from 
the many-particle Schrodinger equation. QHD description 
of unpolarized ultracold fermions was presented in Ref.s 
Gland fl8ll . 

At studying of ultracold electrically polarized Fermi gas, 
many researchers actually prefer using the kinetic equation 
|@1, IToll . Hill , which describes evolution of the distribu- 
tion function /(r, p, t) defined in the six dimensional space 
of coordinate r and momentum p. The kinetic equation al- 
lows us to get a set of hydrodynamics equations for evolu- 
tion of the particle concentration n(r, t) = J /(r, p,t)dp, 
momentum density j(r, t) = J p/(r, p, t)dp, energy and 
higher moments of the distribution function. However the 
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FIG. 2. (Color online) The figure presents the picture of particles 
having electric dipole moment where direction of electric dipole 
moments evaluate in time (rotating or oscillating) and have differ- 
ent directions in different points of space (it might be caused by 
difference of phase at rotation). Arrow by the circle imitates par- 
ticle dipole and it's direction. Round lines around arrow minds 
rotation of dipoles. 



usually used kinetic equation corresponds to the hydrody- 
namic equations, which do not include the equations of po- 
larization evolution, and, so it describe polarization evo- 
lution caused by changing of particles concentration only. 
We believe that such approximation might be enough for 
getting description of collective excitation of the sample 
shape, but it is not enough for bulk waves description in 
the electrically polarized Fermi gas. 

In Ref.s |@], ifioll the kinetic equation was derived via 
the Wigner function with parallel dipoles without dynamics 
of electric dipoles direction. In Ref. Il ill this equation has 
been used for dynamics of the dipolar Fermi gas at finite 
temperature. The functional renormalization group tech- 
nique used in Ref. [19] to get the zerotemperature phase 
diagram of dipolar fermions on a two-dimensional square 
lattice at half filling for the system of dipoles pointing in 
the same direction. 

A. R. R Lima and A. Pelster OOll presented the deriva- 
tion of the continuity and Euler equations for dipoles 
fermions considering the evolution of the one -body den- 
sity matrix l2lll - Ir23ll . They considered the set of equa- 
tions which is equivalents to the corresponding non-linear 
Schrodinger equation and describes dynamics of parallel 
dipoles. In QHD description these equations are only a part 
of the set of QHD equations. Here we present more gen- 
eral set of QHD equations which fully described dynamics 
of electrical polarization in polarized Fermi gas. The Hab- 
bard's model is also very useful lr24ll . lT25ll . but we are not 
going to discuss it here. 

The set of the many-particle QHD equations obtained in 
Ref.s 113I1 - 111511 contains four equation, which are the con- 
tinuity equation, the Euler equation (the momentum bal- 
ance equation), the equation of polarization evolution and 
the equation of polarization current evolution. These equa- 
tions appear along with the equations of field, which are 



VE = — 4-7rVP and V x E = 0, where E is the electric 
field, and P is the density of electric dipole moment, they 
are a pair of the Maxwell's equations. This model give 
us possibility to study the evolution of polarization, and 
it's influence on particle motion. It was found that in lin- 
ear approximation where are two wave solution instead of 
the Bogoliubov's mode existing in the unpolarized Bose- 
Einstein condensate. 

In this paper we develop the many -particles QHD for the 
electrically polarized ultracold fermions. We derive corre- 
sponding the continuity equation, the Euler equation, the 
equation of polarization evolution and the equation of po- 
larization current evolution. Expecting to obtain two wave 
solutions instead of one existing in the unpolarized ultra- 
cold Fermi gas we consider linear approximation of the 
QHD equations. 

This paper is organized as follows. In Sec. II we present 
the derived set of the QHD equations for the electrically 
polarized ultracold fermions. In Sec. Ill we describe the 
method used for the solving of the QHD equations and 
present formula for dispersion dependence of waves in the 
system of the electrically polarized ultracold fermions. In 
Sec. IV we present numerical analysis of the dispersion 
dependence. In Sec. V we present the brief summary of 
our results. 



II. BASIC EQUATIONS 

Here we will briefly present the set of the QHD equa- 
tions for the electrically polarized ultracold Fermi gas. The 
method of theirs derivation described in Ref.s 11411 and 
lfl7ll . It was done for the Bose-Einstein condensate, but 
the method of derivation is the same, this is why we do not 
present details of derivation and present resulting equations 
only. 

The first equation of the QHD equations system is the 
continuity equation 



d t n + d a (nv a ) = 0, 



(1) 



where n is the particles concentration and v a is the velocity 
field. 

The momentum balance equation for the polarized ultra- 
cold fermions has the form 

mn(d t + vVK + —(?>TT 2 )id a ni 
5m 



-—d a An + —d p 

4m 4m \ n 



= -T 2 (d a n)An - 2T 2 (d a d fi n)df 3 n 
+ 3T 2 rad a Ara - 4(37r 2 )§T 2 d Q n3 + P^d a E , (2) 



3 



where 
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(3) 
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In equation (0 we defined a parameter T 2 as pi ). This 
definition differs from the one in Ref.s lfl2h . lll8ll . In left- 
hand side of equation we have three terms proportional to 
H 2 , first of them is the gradient of the Fermi pressure. Other 
two are the quantum Bohm potential, they appear as a result 
of using of the quantum kinematics. The first four terms in 
the right-hand side of equation (0 describe the short range 
interaction in the system of ultracold fermions appearing 
in the third order by the interaction radius. They occur be- 
cause of taking into account of the SRI potential The 
interaction potential XJ^ defines the macroscopic interac- 
tion constant T 2 . The first three terms in the right-hand 
side of equation contain high space derivatives of parti- 
cle concentration and, thus, have something common with 
the quantum Bohm potential. The fourth term arises be- 
cause of the dependence of the interaction on the Fermi 
pressure. Therefore the Fermi pressure gives contribution 
in two terms, kinetic and dynamical. The last term in the 
equation (0 describes force field which affects the dipole 
moments in a unit of volume as the effect of the external 
electrical field and the field produced by other dipoles. The 
last term is written using the self-consistent field approxi- 
mation fi7h . 

We have also field equations 

VE(r,t) = -47rVP(r,t), (4) 

and 

VxE(r,t) = 0. (5) 
Equations (0 and (0 allow us to consider longitudinal 
waves only, i.e. electric field of wave parallel to the di- 
rection of propagation. 

In the case particles does not contain the dipole mo- 
ment, the continuity equation and the momentum balance 
equation form a closed system of equations. When the 
dipole moment is taken into account in a momentum bal- 
ance equation, a new physical value emerges, a polariza- 
tion vector field P a (r, t). This causes system of equations 
to become incomplete. 

We need next equation for investigation of the dispersion 
of the collective excitations is the equation of polarization 
evolution 

d t P a + d R aP = 0, (6) 

R a °(r, t) is the current of polarization. 

The equation © does not contain information about 
the effect of the interaction on the polarization evolution. 
The evolution equation of R al3 (r, t) can be constructed by 
analogy with the above derived equations. Using a self- 
consistent field approximation of the dipole-dipole inter- 
action we obtain an equation for the polarization current 
R ap (r,t) evolution 



d t R af> + 



P a v v 1 



1 fi2 
_|__07 r a07 _ Q AP a 

m 4m 2 
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( dpP a d y n 



8m 2 



+ 



n 



d 1 P a d p n 
n 



a P a P 7 



m n 



(7) 

Here r a ^ is an analog of the tensor of kinetic pressure. 
Assuming to the fact that we consider the system of the ul- 
tracold fermions, and the pressure p described by the Fermi 
pressure (the second term in the left-hand side of equation 
(0). Therefore we suggest following equation of state for 
ther"^ 



r a ^ = K h 2 (37r 2 f^n 2 ^P a 5 fl y5m, 



(8) 



where k is a numerical constant. Functionally r^ 1 de- 



pends on concentration r 



~ n 



2 / 3 and polarization 



~ P a . If we consider a system of particles with par- 
allel electric dipole moments, which direct along z axis, 
we get that P a = p z ^nb za , where p z {) is the constant elec- 
tric dipole moment of single particle. Thus P a changes 
due to change of concentration n only. This approximation 
corresponds to using of the non-linear Schrodinger equa- 
tion. In this case we also get that R a/3 = p^nv^d 20 and 
equation (0 reduces to the Euler equation (0 where in the 
last term we should put P a = p^n5 za . From this compar- 
ison we conclude that k should be equal to 1, but we will 
keep quantity n through the paper to trace contribution of 
the Fermi pressure appearing in the equation of the polar- 
ization current evolution. The last term in the formula (0 
includes both external electrical field and a self-consistent 
field that particle dipoles create. This term contains a nu- 
merical constant a. 

Equation (0 does not contain short range interaction 
because the short range interaction gives no contribution 
in equation ([7]) in the first order by the interaction radius 
(analogously to equation (0), and we do not consider con- 
tribution of higher order. For comparison we can admit that 
in the Euler equation (0 we have accounted the short range 
interaction up to the third order by interaction radius. 

Described equations correspond to the many-particles 
microscopic Schrodinger equation, where for dipole-dipole 
interaction we have used next formula 



H. 



-d a d p 



1 



d«d p 2 . 



Which differs from usually used Hamiltonian of dipole- 
dipole interaction by the term proportional to the Dirac's 
delta function. 

We present and use the set of the QHD equations in the 
form it has been derived. Thus we consider equations of 
field © and Q. In Ref. [16], at studying of electri- 
cally polarized Bose-Einstein condensate, the QHD equa- 
tions was considered along with the whole set of Maxwell's 
equation, what allowed consider transverse waves in polar- 
ized quantum gases. In the result it was shown that electro- 
magnetic waves splits on two branches, and matter waves 
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also contain transverse component that leads to anisotropy 
of the spectrum of collective excitations, whereas using of 
equations I© and Q allows to consider longitudinal waves 
only what gives no anisotropy. 



III. ELEMENTARY EXCITATIONS IN THE 
ELECTRICALLY POLARIZED FERMI GAS 

We can analyze the linear dynamics of the collective ex- 
citations in the polarized BEC using the QHD equations 
CQ- ©• in the beginning we consider the system is placed 
in the external electrical field parallel to direction of wave 
propagation E = E e x . Below we will consider influence 
of the equilibrium external electric field E 0z = E 0z e z , so 
we will have E = E 0x e x + E Qz e z . The values of concen- 
tration n and polarization P = kE for the system in the 
equilibrium state are constant and uniform and its velocity 
field v a (r, t) and tensor R a/3 (r, t) values are zero. 

We consider the small perturbation of equilibrium state 
like 



n = n + 5n, 



+ v a , E a 



Eg + 5E a 



pa = pa + Spa^R"? = + 5R aP . 



(9) 



Substituting these relations into system of equations £[]l- 
$Q and neglecting nonlinear terms, we obtain the set of 
linear homogeneous equations in partial derivatives with 
constant coefficients. Passing to the following representa- 
tion for small perturbations 5f in the form of plane wave 
propageting in direction of x axis 

5f = f(uj, k)exp(—iut + ikx) 

yields the homogeneous system of algebraic equations. 
The electric field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system 
of equations in terms of the electric field, we come to 
dispersion equation. Equation dU gives us the dispersion 
equation and equation ([5]l gives us additional condition on 
8E y = 0, 5E Z = 0, so all quantities we express via 5E X . 
Thus, in presented model we have deal with longitudinal 
wave as self-consistent electric field in the wave parallel 
to the direction of wave propagation. Solving this equation 
we get despersion dependesies. Here we discuss dispersion 
equation emerging from evolution of 8n, 5v x , 5P X , 5R XX 
and 5E X . 

The dispersion characteristic for collective excitations in 
the electrically polarized Fermi gas can be expressed in the 
form of 




+ (1 + ^)^(3^)2/3^/3^ 



| 3T 2 re fc 4 + AnaP^k 2 + 32 (37r 2 n 
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47raP 2 fc 2 + 32 /3 no fc2 ^ 
3 m 



mn 



+ 32ttk(3tt 



2\2/3. 



n 2 



k^p 2 



1/2 N 



15m 2 mn Q 



(10) 



The last term under the square root, which is the last term in 
formula ([Tol l, is positive. Other term under the square root 
is the square of sum of terms, so we have the sum of two 
positive quantities under the square root. This means that 
LO 2 has no imaginary part for considered system of particles 
and we can conclude that there are no linear instabilities in 
this case. 

For numerical analyzes of obtained formula ( fTOt we 
introduce dimensionless parameters Q, = mio/fon^ 3 ), 

7 = mn T 2 /Ji 2 , d = mP^ /(h 2 n^ 3 ), and £ = k/n^ 3 . 
Let's rewrite formula (TTOt in dimensionless variables 



Q 2 



e 



1 ^2 3 , 2 1/1 Hi 

+ 2^ + 2 3 + 5 



(3^) 



2\2/3 



+W + y (3vr 2 ) 2 / 3 7 ± ^((3vr 2 ) 2 / 3 Q - f 



_ 4w(i+ £f(3^ ) 2/3 7 + 3 7 ^ 
3 



+ ^( 3 ^ 2 ) 2 /^ 2 
15 



1/2 N 



(ID 



The first term in this formula appears from the quantum 
Bohm potential, the second and nine terms exist due to SRI 
which we account up to the third order by interaction ra- 
dius. The third and sixth terms are caused by the Fermi 
pressure which is a part of equations (O and (O via kinetic 
terms. The fourth and seventh terms present contribution 
of equilibrium polarization. The fifth and eighth terms ap- 
pear as consequence of contribution of the Fermi pressure 
in the SRI. The last term in formula ( fTTT) appears due to si- 
multaneous account of the equilibrium polarization and the 
Fermi pressure. 

First of all we should consider the limit of a small polar- 
ization contribution, to compare it with dispersion of unpo- 
larized quantum Fermi gases. Making series of the square 
root in formula ([Til on small parameter d up to the linear 
terms we get 



n 2 



i 



(37T 2 ) 



2\2/3 



1 +32- )r + ( \ + 3 1 



5 



8-7TK 



(3TT 2 ) 2/3 d 



15 



(37T 2 ) 2 /3 i - I + f (3tT 2 ) 2 / 3 7 + 3^2 



(12) 



for the solution with the plus infront of the square root, and 
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(37r 2 ) 2 / 3 d 
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(13) 

for the solution with the minus infront of the square root. 

For comparison we present here the dispersion of plane 
wave in unpolarized ultracold Fermi gas in introduced 
above dimensionless variables 

n 2 = i(3vr 2 ) 2 / 3 (l + 32 7 )e 2 + Q + 3 7 )C 4 , (14) 

where following paper 112I1 we include contribution of SRI 
between Fermi particles up to the third order by the inter- 
action radius. 

We can see that solution (fT2l differs from (fT4l) by one 
term, the last one proportional to d. In the case of par- 
ticles having no electric dipole moment d = we have 
that the solution ( TTTT ) with plus infront of the square root 
describes well-known matter wave in ultracold Fermi gas, 
and formula ( fT3l) presents dispersion dependence of new 
wave appearing due to polarization dynamics. 

Introducing the average distance between particles I = 
1 / ^/no we can rewrite definition of £ in the following form 
£ = 2nl/\, where A is the wavelength A = 2ir/k. As 
wavelength must be at least more than two distances be- 
tween particles we have that £ < 1 or, in more realistic 
cases, we should consider £ ~0. 1-0.001. Thus we can 
admit that £ 2 on several orders smaller than one. To get 
7 = 1 we need | T 2 |> 10~ 20 ^/JIg(fim) 5 5 /c. For equi- 
librium particle concentration equal to 10 14 cm -3 , particle 
mass m = 2 ■ 10~ 23 g and electric dipole moment of par- 
ticle p = 1 Debuy we have d = 0.8. More interesting for 
future experiment case when molecule has electric dipole 
moment equal to 100 Debuy, in this case d ~ 10 4 . For- 
mula (fTTT l presents two solutions, for the case when parti- 
cles have large electric dipole moment for solution with the 
sign plus in front of square root we get 



(15) 



or in the other terms u = h^/n^V Airadk j m, this solu- 
tion appears at neglecting all terms in comparison with two 
terms proportional to ad. Let's consider the second solu- 
tion presented by formula (fTTTl with sign minus in front of 
the square root. If we left two terms only, which propor- 
tional to ad we get $7 = 0. Thus, we have that two large 
terms reduce each other, so to get the solution we should 



consider other terms. Making expansion of the square root 
in linear on 1/d approximation we get dispersion of the 
second mode 



tt 2 



e (I(3vr 2 ) 2 / 3 (l + :S2, 



+ 3 7 



2 K 
15a 



(3vr 2 ) 2/3 ]£ 



(16) 



From comparison of formulas (fT4l and (fT6l l we can see that 
in the limit of a large equilibrium polarization dispersion 
dependence of the "second" matter wave similar to the dis- 
persion of matter wave in unpolarized ultracold Fermi gas 
(fT4l . with only difference - existence of the last term in 
formula (fTol l. 

Following Ref . ifl^l we can admit that account of contri- 
bution of the transverse electric field in dynamics of matter 
waves in electrically polarized ultracold Fermi gas should 
lead to replacement of d by dcos 2 9. Thus, we can write 
that in general case we get 



O = V^irad • | cos0| • £ (17) 

instead of (fl~5T > for dispersion of the collective excitations 
in electrically polarized ultracold Fermi gas of molecules 
with large electric dipole moment. We can also admit that 
made replacement has no influence on solution ( TToT) . 

Here we have considered dispersion dependence appear- 
ing from dynamics of Sn, 5v x , 5P X , 5R XX and 5E X . Quan- 
tities 8P y , 5R yx and 5P Z , 5R ZX give us two pair of inde- 
pendent sets of equations. Both of them give same disper- 
sion relation: 



h 2 k 4 



CO 



+ K(37T 2 n 



,2/3 



tr 



:k 2 



Am 2 5m 2 
Taking into account E 0z ^ (or E 0y ^ 0) gives no 
influence on evolution of 5n, 5v x , 5P X , 5R XX and 8E X . 
However it gives contribution in evolution 5P Z , 5R ZX . 5n 
gives influence on 5P Z , SR ZX , but in this case we have not 
closed set of equations. 

Now, we have described main properties of dispersion of 
collective excitation and we can pass to conclusions. 



IV. CONCLUSION 

We have developed the method of the many-particle 
QHD for the system of the ultracold electrically polarized 
Fermi particles. The set of the QHD equations consists of 
four material equations: the continuity equation (the equa- 
tion of balance of the particles number), the Euler equation 
(the momentum balance equation), the equation of polar- 
ization evolution, and the equation of polarization current 
evolution, and two equations of field, which is the part of 
the set of the Maxwell's equation. Due to the fact of using 
of the mentioned equations of field we can admit that we 
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FIG. 3. (Color online) The figure describes the dispersion depen- 
dence of the collective excitation in the polarized ultracold Fermi 
gas presented by formula l |17| | at d = 100 and a = 1. 

have deal with the longitudinal waves, i.e. direction of the 
electric field in the wave parallel to the direction of wave 
propagation. 

Using the QHD equations we have studied dispersion of 
the plane collective excitation in the three dimensional ul- 
tracold electrically polarized Fermi gas. We have found 
that there are two matter waves, whose dispersion reveals 
no anisotropy. These two waves appear in the electrically 
polarized Fermi gas instead of the one matter wave exist- 
ing in the unpolarized Fermi gas. Using our previous re- 
sults we generalized the obtained dispersion dependencies 
to get anisotropic spectrum. 

ACKNOWLEDGMENTS 

The author thanks Professor L. S. Kuz'menkov for fruit- 
ful discussions. 



[1] M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, 

arXiv: 1207. 1914. 
[2] Patrick Koberle, Holger Cartarius, Tomaz Fabcic, Jorg Main 

and Gunter Wunner, New Journal of Physics 11, 023017 

(2009). 

[3] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein and T. 

Pfau, Rep. Prog. Phys. 72, 126401 (2009). 
[4] Lincoln D. Carr, David DeMille, Roman V. Krems and Jun 

Ye, New Journal of Physics 11, 055049 (2009). 
[5] K.-K. Ni, S. Ospelkaus, D. J. Nesbitt, J. Ye and D. S. Jin, 

Phys. Chem. Chem. Phys. 11, 9626 (2009). 
[6] S. Yi and L. You, Phys. Rev. A, 61, 041604(R) (2000). 
[7] K. Goral, K. Rzazewski, and T. Pfau, Phys. Rev. A 61, 

051601(R) (2000). 
[8] L. Santos, GV. Shlyapnikov, P. Zoller, and M. Lewenstein, 

Phys. Rev. Lett. 85, 1791 (2000). 
[9] T. Sogo, L. He, T. Miyakawa, S. Yi, H. Lu, and H. Pu, New 

J. Phys. 11, 055017 (2009). 
[10] J.-N. Zhang, and S. Yi, Phys. Rev. A 80, 053614 (2009). 
[11] J.-N. Zhang, R.-Z. Qiu, L. He, and S. Yi, Phys. Rev. A 83, 

053628 (2011). 

[12] P. A. Andreev, L. S. Kuz'menkov, Phys. Rev. A 78, 053624 
(2008). 

[13] P. A. Andreev and L. S. Kuz'menkov, arXiv: 1106.0822. 
[14] P. A. Andreev and L. S. Kuz'menkov, arXiv: 1201.2440. 
[15] P. A. Andreev, Russian Physics Journal 54, 1360 (2012). 
[16] P. A. Andreev and L. S. Kuz'menkov, arXiv: 1208. 1000. 
[17] P. A. Andreev, L. S. Kuz'menkov, M. I. Trukhanova, Phys. 

Rev. B 84, 245401 (2008). 
[18] K. V. Zezyulin, P. A. Andreev, and L. S. Kuzmenkov, 

arXiv:1205.1161. 
[19] S. G. Bhongale, L. Mathey, Shan- Wen Tsai, Charles W. 

Clark, and Erhai Zhao, Phys. Rev. Lett. 108, 145301 (2012). 
[20] A. R. P. Lima, A. Pelster, Phys. Rev. A 81, 063629 (2010). 
[21] G. M. Bruun and C.W. Clark, Phys. Rev. Lett. 83, 5415 

(1999). 

[22] M. Amoruso, I. Meccoli, A. Minguzzi, and M. Tosi, Eur. 

Phys. J. D 7, 441 (1999). 
[23] P. Ring and P. Schuck, The Nuclear Many-Body Problem 

(Springer, Berlin, 2004). 
[24] Zhen-Kai Lu, G. V. Shlyapnikov, Phys. Rev. A 85, 023614 

(2012). 

[25] Anne-Louise Gadsbolle, and G. M. Bruun, Phys. Rev. A 85, 
021604 (2012). 



* andreevpa@physics.msu.ru 



